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Having in mind its future extension for theoretical investigations related to charmed nuclei, we 
develop a relativistic formalism for the nonmesonic weak decay of single-A hypernuclei in the frame¬ 
work of the independent-particle shell model and with the dynamics represented by the (rr, K) 
one-meson-exchange model. Numerical results for the one-nucleon-induced transition rates of 
are presented and compared with those obtained in the analogous nonrelativistic calculation. There 
is satisfactory agreement between the two approaches, and the most noteworthy difference is that 
the ratio r^/Tp is appreciably higher and closer to the experimental value in the relativistic calcu¬ 
lation. Large discrepancies between ours and previous relativistic calculations are found, for which 
we do not encounter any fully satisfactory explanation. The most recent experimental data is well 
reproduced by our results. In summary, we have achieved our purpose to develop a reliable model 
for the relativistic calculation of the nonmesonic weak decay of A-hypernuclei, which can now be 
extended to evaluate similar processes in charmed nuclei. 

PACS numbers: 21.80.-|-a, 13.75.Ev, 21.60.-n 


I. INTRODUCTION 

Investigations of exotic nuclear properties, such as 
large isospin (manifest in the so called neutron-rich iso¬ 
topes), or nontrivial values of flavor quantum numbers 
(strangeness, charm or beauty), are of continuous inter¬ 
est. The best known nuclei within the last category are 
those where a A-hyperon, with strangeness S = — 1, re¬ 
places one of the nucleons, giving to the composed system 
some quite unusual properties. Such nuclei are referred 
to as A-hypernuclei — for recent reviews, see Refs. [l|,[3- 

One of the most remarkable properties of A-hyper- 
nuclei is the occurrence of the nonmesonic weak decay 
(NMWD), induced by the elementary process 

A + N ^ n + N, (1) 

with N = p (proton) or n (neutron). This is the main 
decay channel for medium- and heavy-weight hypernu¬ 
clei — Refs. 3 and 0 provide, respectively, reviews 
on recent theoretical and experimental developments in 
the study of hypernuclear decay; for earlier comprehen¬ 
sive reviews on theory see Refs. [ShTI and on experiment 
Refs. [1,1^. NMWD can only take place within the nu¬ 
clear environment and is a unique opportunity offered by 
nature to access the strangeness-changing interaction be¬ 
tween baryons. Its mean lifetime has been measured in 
several A-hypernuclei and found to be of the same order 
of magnitude as the full mean lifetime of A in free space, 
ta = (2.632 ± 0.020) x 10-i° s [I^. 

The NMWD dynamics is frequently handled by one- 
meson-exchange (OME) models. Such models are moti¬ 
vated by the fact that the NN interaction at long dis¬ 
tance is due to the one-pion-exchange, but with the differ¬ 
ence that in NMWD the exchange processes occur with 
one strong and one weak vertex and can include other 


mesons in addition to the tt, like the pseudoscalar {K, rf) 
and vector (p, w, iT*) mesons [lll - l^ . The coupling con¬ 
stants at the strong vertices can be taken from different 
OME models for the NN interaction, while those at the 
weak vertices can be extracted from free A decay data 
and making use of soft meson theorems and SU{6)w 
symmetry [111 . Il^ . A recent study conducted within a 
nonrelativistic framework [2^ indicates that tt and K 
exchanges give the main contributions to the NMWD of 
s-shell hypernuclei. 

Instead of implanting a A in a nucleus one could also 
imagine to implant a charmed baryon, like e.g. a A)!", 
in view of the similarity between the quark structures 
of the strange and charmed hyperons, namely A{uds) 
and N^{udc). Such a possibility was in fact conjec¬ 
tured 40 years ago [2^ and several authors in the suc¬ 
ceeding decades have found, using different models for 
the interactions between nucleons and charmed hyper¬ 
ons, that such hypothetical exotic nuclei (including even 
bottom nuclei) could actually form a rich spectrum of 
bound states over a wide range of atomic numbers [2^ 
Id^l . Like A-hypernuclei, A+-hypernuclei may also decay 
via a NMWD process. One example is [4l| 

A+ -I- n A -I-p, (2) 

which can be induced by the exchange of a tt, p or AT 
meson. Another possibility is 

At + N ^p + N, (3) 

induced by the exchange of a D meson. Experimentally, 
the literature only reports, inconclusively, the formation 
of three A)!'-hypernuclei, observed in a series of emulsion 
experiments [dj, |4^. But this situation can change in 
a few years, with the starting of operation of the FAIR 
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facility in Germany and the Hadron Facility at JPARC 
in Japan. 

There are, however, important differences between 
NMWD in A-hypernuclei and A+-hypernuclei. A first 
difference comes from the mean lifetimes of the two hy- 
perons: r^+ ^ 10“^ r^. While the mean lifetime of the 
NMWD (IJ) is of the same order of magnitude of the 
full mean lifetime of A in free space, no theoretical esti¬ 
mate has been made for the decays @ and (1^ . In ad¬ 
dition, while the free-space decay of A is dominated by 
the pionic channels A —^ p7r“ and A —^ mr^, with other 
decay channels contributing a thousand times less, A+ 
decays in two semileptonic and numerous hadronic chan¬ 
nels with S = — 1 final states, having branching ratios 
of a few percent each. Also, decays into channels with 
S = 0 and S = -1-1 are Cabibbo-suppressed by factors of 
the order of 10“^-10“^ M- A second very important 
difference concerns the energy liberated in the decays, 
which is of the order of the mass difference A of the 
particles involved in the weak vertex: for the decay @, 
A = — M\ = 1170.9 MeV, and for the decay ([3]), 

A = M^+ — M^v = 1348.2 MeV, which should be com¬ 
pared to A = Ma — Mat = 177.3 MeV for the decay ([T]). 
One consequence of such large energy releases is that non- 
relativistic approaches, like those of Refs. [TiI - It^ . be¬ 
come inapplicable for the evaluation of NMWD transition 
matrix elements in charmed hypernuclei. In addition, a 
large energy release also implies that nuclear recoil can¬ 
not be neglected in the calculation of decay rates, par¬ 
ticularly for light-weight nuclei. On the other hand, the 
interactions of the fast outgoing nucleons and/or hyper- 
ons with the residual nuclear system are expected to play 
a minor role. 

In the present paper we develop a relativistic formalism 
for NMWD of hypernuclei within an independent-particle 
shell model (IPSM), and discuss the inclusion of recoil. 
Although the use of a relativistic model for the study of 
the structure of hypernuclei dates back to the late 1970’s 
with Brockmann and Weise [d^, so far little is known 
about the impact of a relativistic approach in the evalu¬ 
ation of NMWD rates. The first studies started 25 years 
ago with Ramos et al. [4^,|4^. These authors used single¬ 
particle bound-state wave functions obtained by solving 
the Dirac equation with static Lorentz-scalar and -vector 
Woods-Saxon potentials, and transition matrix elements 
calculated with a (tt, K) OME model. More recently, a 
similar approach was used by Conti et al. [U13 > where 
the nuclear structure was described by a finite-nucleus, 
relativistic, mesonic-mean-field model. An interesting 
feature of these studies is that the reported numerical 
results for the decay rates differ considerably from those 
obtained with nonrelativistic approaches — this is true, 
e.g., for the hypernucleus, as we show in Table Hill 
Such differences are larger than one would expect given 
the moderate energies involved in the decay process. It is 
also important to notice that these predictions strongly 
contradict the experimental data. 

Our aim in the present paper is to set up a relativis¬ 


tic formalism for NMWD with the perspective of future 
applications to charmed hypernuclei. In other instances 
involving nuclear structure calculations at low and inter¬ 
mediate energies, it is often more convenient and sim¬ 
pler to use a relativistic approach than a nonrelativistic 
one [43; this seems to be also the case for NMWD — 
Ref. [SOjl presents a very complete review on relativistic 
approaches for the study of nuclear structure. Although 
our approach for the NMWD of hypernuclei shares sim¬ 
ilarities with the formalism of Refs, there are 

noteworthy differences: 

1. Our final expressions for the decay rates do not in¬ 
volve angular momentum projection quantum num¬ 
bers, since they have been summed over in closed 
form using the Racah algebra, which simplifies the 
numerical calculation; 

2. Spectroscopic factors are evaluated in the sec¬ 
ond quantized formalism, as done for instance in 
Ref. [^, without recurring to the technique of co¬ 
efficients of fractional parentage (c.f.p.’s), which is 
the standard antisymmetrization procedure in the 
first quantization framework, see e.g. Ref. [s^ : 

3. We discuss the inclusion of recoil. 

The predictions of our formalism are compared with 
available data [53 - 3 ^r the NMWD rates of the 
hypernucleus. In addition, we make a detailed compar¬ 
ison with results obtained in nonrelativistic approaches 
that include the same ingredients (like short-range cor¬ 
relations and OME model); such a comparison between 
the outcomes of analogous relativistic and nonrelativistic 
approaches had not been done so far. 

Our formalism is explained in Section [IT] starting from 
the simplest scenario, corresponding to hypernuclei with 
closed-shell cores and ignoring recoil, in Subsection III Al 
This part is done in a strictly relativistic manner, while 
the next two steps are performed in analogy to nonrela¬ 
tivistic calculations: first, in Subsection III B[ we general¬ 
ize the formulation to hypernuclei with open-shell cores; 
and secondly, in Subsection IB Cl the recoil effect is dis¬ 
cussed. Subsequently, in Section imi our numerical re¬ 
sults for the decay rates of a^C are presented and com¬ 
pared to those of nonrelativistic calculations using a sim¬ 
ilar model • They are also compared with those 

of previous relativistic calculations and confronted with 
the experimental data, and a few conclusions are drawn. 
Einally, in Section lYl a general summary is given. In 
Appendices lAl - ICl some details of the calculation are pre¬ 
sented. 


II. RELATIVISTIC DECAY RATE 

To derive the NMWD rate we start from the Fermi 
Golden Rule. For a hypernucleus in its ground state with 
spin Ji and total rest energy Ej decaying into (i) two free 
nucleons, with asymptotic kinetic energies (Ti,T 2 ), spin 
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projections (si,S 2 ), and isospin projections (^ 1 ,^ 2 ) and 
(ii) the residual (^4 — 2)-system, with spin Jp, total rest 
energy Ep, and kinetic energy of recoil T/j, reads 


r 


nm 


51-52^1^2 


dpi dp 2 
(27r)3 (27r)3 


S{Ei 


Ep-£) 


coupled to an (A — 1 ) nuclear core of spin Jc and total 
rest energy Ec, i-e., |J/) = \{JcjA)Ji), having energy 
Ej = Ec + Eja + Ma; (ii) the nucleon N inducing the 
decay is in the single-particle state jat (j = nlj)] (iii) 
the final residual nuclear states have the form \Jp) = 
\iJcjN^)JF) with energy Ep = Ec - - Mat; (iv) the 

liberated kinetic energy is 


X \A4{piP2SlS2tlt2JFMp, , 


(^) TP+T2+T1 — Ei — Ep — 2 Mn = Ajat — A-l-EjA+Ejiv (d) 


where £ = 2MAr - Tr - T 2 - Ti, M = {1 - Pi2)MlV2 
is the antisymmetrized and normalized relativistic ma¬ 
trix element that is specified below, Mat is the nucleon 
mass, and pi = \/Ef — and Ei = Ti + Mm are the 
asymptotic momenta and total energies of the outgoing 
particles {i = 1,2). We use unitary, as opposed to co¬ 
variant, normalization for the momentum eigenspinors; 
for details see Section 2.2 of Ref. [i^- We average over 
the spin projections M/ of the initial hypernucleus and 
sum over the final spin projections Mp. For the nuclear 
structure framework, the IPSM is used, while the dynam¬ 
ics is described by an OME potential containing always 
one weak vertex W and and one strong vertex S, as illus¬ 
trated in Figure [TJ In the IPSM it is assumed that: (i) 



FIG. 1. Diagrammatic representation of the hypernu- 
clear nonmesonic weak decay from the initial state \Ji) = 
\{JcjA)Ji) to the final state |Jf) = \{JcjM^)JF) while two 
nucleons with momenta pi and P 2 are emitted into the contin¬ 
uum. S and W are the strong and the weak vertices, respec¬ 
tively, and M is a nonstrange meson. For a strange meson, 
the natures of the two vertices should be interchanged. 

the initial hypernuclear state can be approximated as a 
A-hyperon in the single-particle state j’a = lsi /2 weakly 


where A = Ma — Mat = 177 MeV, and the e’s are single¬ 
particle energies. 

A. Hypernuclei with doubly-closed shell cores and 
without recoil 

Taking the simplest possible case in Eq. o, we will 
start with hypernuclei whose cores contain only doubly- 
closed subshells, as, for instance, aHc, )J0 , and we 
will omit the recoil effect. Thus, Jc = 0, Jj = j'a, Mi = 
WA, Jf = JN-, Mp = mM, and the transition amplitude 
M. is just the two-body T-matrix for the direct OME 
process. When a pseudoscalar coupling is considered for 
the strong vertex, one has for the pio n plus kaon meson 
exchange: M = M'" + , see Eq. (3)] and [13, 

Eqs.(7),(45)], with 

{piP2SlS2tlt2jNmMjKmR) = 

J dxdy'(/ip,sj(x)[yl’"(ti,t2) - 12 ) 15 ] 

X ^jArriA (x) A’^dx - y\)^p^s 2 {y)l5'i’jjv7nM (y): (6) 

and 

{piP2SiS2tit2jNmNjAmA) = 

J dx dy (x)75 4'jAmA (x) A*'(|x - y|) 

X'0p2S2(y)[-^^(^1.^2) - B^{ti,t2)l5]'i’jNmr,{y), (7) 

where we are using the following definitions 
{ti,t2) = CpTn^g-^MN-M {ti,t2), 

A^{ti,t2) = G pm^gKAN A^ iti,t2), 

B^{ti,t2) = Gpmlg^MNB^{ti,t2), 

B^{ti,t2) = GpmlgKANB^ {ti,t2), ( 8 ) 

with A'^{ti,t 2 ) = AI, A^{ti,t 2 ) = {lAi + KAo), 
B^{ti,t 2 ) = BG B^ih,t 2 ) = (IBi + KBo) and 

where I = (ti|Ti|tA = -^) • {t 2 \T 2 \tN) and K = 
(ti|li|tA = — 5 ) (^ 2 |l 2 |iAf) are, respectively, the isovector 
and isoscalar isospin factors. Here, Gpml. = 2.21x 10“^, 
with Gp being the Fermi weak constant and the pion 
mass, while g-^NN = 13.3 and gxAN = —14.1 are the 
strong vertex couplings [ 6 l[- The pion parity-violating 
(PV) and parity-conserving (PC) weak coupling con¬ 
stants are adjusted to the free A-decay giving, respec¬ 
tively, A = 1.05 and B = —7.15, while the kaon weak 
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couplings 


riPV nPV 

= A^= ^ 


Bn = 




^ B, = 


2 


2 • ~ i\ 7 ~ i 2 ’ 

with = 0.76, = -18.9, OP' = 2.09, and 

Dp^ = 6.63, have been estimated theoretically lllll. 

The propagator A^(|x—y|) (M = n,K) reads |4^ . l4^ : 

1. For (go)^ < mh, 


A"(|x-y|) = - 


exp(-- 9o |x - y|) 


47 r|x-y| 


+ A^(x-y), 


( 10 ) 


where qo is the energy carried by the exchanged 
meson, and 


,M/ 


-y|) = - 


exp(-\/A^-g^|x-y|) 


Am- m 


2 

M 


Sttv'A^ - 


Qq 


47 r|x-y| 
exp{-^jAlj - ql |x-y|), ( 11 ) 


is the finite-size correction when dipole form fac¬ 
tors with cut-off parameters Am {M = tt,K) are 
attached at each vertex; 

2. For mh < {qof < Aj^, 


_aM/|„ ..in exp(iv^gg-TO^|x-y|) 

47 r|x-y| 

Therefore the propagator is now complex [13, 
and can have oscillatory behavior, with real and 
imaginary parts given by 


5RA^(|x-y|) = - 

|x-y|) = - 


5A^( 


cosjy'qo - m;^|x-y|) 
47 r|x-y| 

sin(\/g^ -TO^|x-y|) 


47 r|x- y| 
in this region of transferred energy. 


(13) 

(14) 


The state of each ejected nucleon, with asymptotic mo¬ 
mentum p and spin projection s, will be approximated 
by a Dirac plane wave, which is expanded in spherical 
partial-waves as follows [13, Appendix D]: 


V’ps(r) = y^/pslKm}*-ippKm(r), 

Km 

{ps\Km)* = AW'^{l^i^s\jm)Yi*pp), 


-iHpK.mi'r) 


(15) 

(16) 

(17) 


where the radial partial-waves are, in unitary norma¬ 
lization. 


fpPr) = \l ji,.ipr) 


(18) 


E — M^v 

gpA'r) = -sgn[K)\l —— jlAP^)7 (19) 


2E 

with K = ±1, ±2 ,..., = |k| - 1 / 2 , 

K for K > 0 


L = 


— K—1 for K < 0 ’ 


( 20 ) 


and Ik = 1-k- (To change to covariant normalization, 
used in Refs. [46l - l4^ . make the replacement —>• 

•\/2M)v in Eqs. m and (HU) and insert the factor 
/{E 1 E 2 ) in Eq.dH).) The angular part is written, in 
standard notation, as 


$«m(r) = ^(;Ai:^s|jm)TiM(r)x^ 


( 21 ) 


and the expansion coefficients (ps|Km)* fulfill the follow¬ 
ing relations 

/ dp{ps\Km)*{ps\K'm') = (ATrfSKK'Smm', ( 22 ) 

S 

2pSjj>^^ j dp {ps\Krn)* {ps\k' m) ■ ■ ■ = 


(47r)^5« 


d cos 9 ■ 


(23) 


where we are using the notation j = y/2j -|- 1. The first 
of these relations can be easily verified, while the second 
one is shown in Appendix The bound-state, single¬ 
particle, wave functions read 


Kjn (l*) — 


1 


Bk{t Km (1*) 


(24) 


r \^-zGK(r)(f)$_«,m(f) 

-A'i'Kmir) 

As explained in Appendix m they are evaluated as in 
Ref. 0 Eq.(16)]. 

To simplify the presentation of formulas in the analyt¬ 
ical development of Eq. 0, the intermediate steps will 
be exhibited only for , which is rewritten as 

A4^(piP2SiS2Di2jArm7v/AmA) 

= (Pl5l|KlWl)(P2S2|zt2TO2) 

Kimi 

K2m2 

X {piKimitiP2K2m2t2\^^\KA'mAKNmN), (25) 

where we are using the following compact notation 
(piKimitiP2K2Tn2t2\^^ iKAmAKNmN) 

= J dxdy'ipp^Kimii^)[A^{ti,t2) - B^{ti,t2p5] 

(xjA’^dx- y|)V’p2K2m2(y)75 4'jNmiv(y)- (26) 
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Introducing these expansions in Eq. 0 gives rise to aux¬ 
iliary quantities such as 

S'"{piti,P2t2) = ^ I dpidp2S{Ajj^-Ti-T2-Tii) 

SlS2 

X |A^’"(piP2SlS2tlt2jJVTOJvjATOA)P, (27) 

in which we evaluate all the summations over angular 
momentum projection quantum numbers and angular in¬ 
tegrations. Neglecting recoil, i.e., setting Tr = 0, we can 
use (1^ for both outgoing particles getting 

S^Plh,p2t2) = {^7T)^ Y, E - ^ 2 ) 

mArriN Kim-i 
K277T.2 

X \{piKimitiP2K2m2t2\A'^\KAmAKi\[mN)\'^ ■ (28) 

Now we perform the angular momentum couplings J = 

I 


JA+JAT and J' = ji +j 2 - As is rotationally invariant, 
it turns out that J = J', which leads to 

{piKlTmtiP2K2m2t2\A'^\nAmAKN'mN) 

= y^(piKitiP 2 A^ 2 ^ 2 ./|A^|KAAjvJ) 

JM 

x(jimij2W2| JM) ijAmAjNmN\JM), (29) 

and 

S^{piti,p 2 t 2 ) = (4^)^ ^ - Ti - T 2 ) 

K 1 K 2 J 

X \{piKitiP2K2t2J\A^\KAKNJ)f , (30) 

where the coupled matrix element of the pion propagator 
is explicitly given by 


{piKitiP2K2t2J\A^\KAKArJ) =-^ j dy | [ A’" (ti, f2) (Vpi «i (x)^ 4 '«A (x) - (x)'^d'*^ (x)) 

-b S^(ti,t 2 )( Vpi«;i(x)-‘'4'«a(x) +'‘-V'pi«i(x)'^4'«a(x)) ] A’^(lx-yl) 

^ (Vp2K2(y)'^^«iv(y)+ '‘'V'p2«2(y)^^«iv(y))}, ,,> (3i) 

{J1J2\JAJN-,J) 


with the above mentioned angular momentum couplings 
indicated in the last index. 

At this point it is convenient to perform the tensor ex¬ 
pansion of the propagators A'^(|x — y|) in the way done 
by de-Shalit and Talmi [H, Sec. 21] for two-body inter¬ 
actions, i.e., 


A-(|x - y|) = ^ Alix, y)[YLi±)L ■ Y{y)L] (32) 

L 


where 


^Lix,y) = 27r 


A’^dx - y|)PL(cos6'2:y)d(cos6>2,y) (33) 


and 


(KlK2d|[YL(x) • YL(y)]|KA/Cjvd) = (-)-?2+M-SJ 

jIl} i^^\\YLMi^2\\YL\\nM). (34) 

It is then easy to demonstrate that 

{KiPitiK2P2t2J\A''\KAKNJ) = + 

L 

X I ? i \ {KlPltlK2P2t2\Al\KAKN) , (35) 

[JN JA ^ ) 


where 


{KiPitiK2P2t2\Al\KAK.N) 

= j xydxdy[B^{ti,t 2 )BYAi.^P^) 

-iA^{ti,t2)AYS^Pi)\^l{x,y)CYAyP-2) 

with 

ALa(^7') = [/p«(»')A«a(?’) - 5p«(?')G«A(r)] 

x(k||Yl||ka), 

BYi^P) = [/p«( 0 G'«a(?’) +£'p«(?’)^«a( 0 ] 
x(-K||yL||KA), 

CL^i^P) = [/p«( 0 G'«a ,(0 +£'p«(?’)^«iv( 0 ] 

x(-K||lL||KAr). 

The reduced matrix elements 

(k||Yi,||k') = (47r)-^/2(-)^-iAjj'L 

f J L f\ 1 + (-)W'+L 

l-i 0 ij 2 


{-k\\Yl\\k’) = (47r) 

( j L f\ 1 + (-)W'+L 

V -1 0 5 J 2 


(36) 


(37) 


(38) 


(39) 


fulfill the symmetry relations (A||yL||K') = (A'||yi||^); 
and (kIITlII - k') = (-k||Yl||k'). 
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The K meson is incorporated through the substitution 

Al = AJ + A^ in (1551) . with 

(«;iPitiK 2 P 2 t 2 | Af IkaKat) 

xydxdy Bj^^^^{xpi)A^{x, y) 

x[B^(ti,t2)C^^^^{yp2) - iA^ih,t2)Df;^^^{yp2)], (40) 
where 

DL^{rp) = [/p«(r)F„„(r) - gp^(r)G^„(r)] 

x(«:||rL||«;jv). (41) 


and get 

= - E ^ ^(A,„ - Ti - T2) 

K,iK2J 

X p(Ti,r2) |(KiPitiK2P2i2d|A|KAKArJ)|^ , (46) 

where 

p(ri,T 2 ) = (Mjv + Ti)A/ri(2Mjv + Ti) 

X (M^ + T2)^yT2{2MN + T2). (47) 


Clearly the above substitution must be accompanied by 
the replacement A'^ —>■ A = A^ + A^ in (1351) , giving 

(AtiPitiK2P2l2d|A|AtAK7vd) = y^(-)J^+JA + J 

L 

X \ "1} {KiPitiK2P2tl\AL\KKKN)- (42) 

[JN JA ^ ) 

Finally, from (|^. 

Fjv = - E 4" / Pi^PiP2dp2 6iAj^ - Ti - T 2 ) 

KIK 2 J 

I _ . 2 

X |(KiPitiK2P2i2'/|A|A{AKArJ)| , (43) 

where 

\K 1 p 1 t 1 K 2 P 2 t 2 J) = ■^( \n 1 P 1 t 1 K 2 P 2 t 2 J) 

- \K2P2t2KlPltiJ)) (44) 

stand for the antisymmetrized and normalized two- 
particle wave functions with the isospins included. The 
isospin factors for the direct and exchange terms of the 
matrix-element in Eq. (l4^ are listed in Table HI 


TABLE I. Isospin factors the for direct (D) and exchange (E) 
terms of the matrix-element in Eq. (1151) . 




I 



K 



n 


p 

n 


P 

D 

1 


-1 

1 


1 

E 

1 


2 

1 


0 


It is worth noting that the matrix elements 
{KipitiK, 2 P 2 t 2 J\A\K,AK,NJ) are in general complex, as 
seen from Eqs. dsnu and (gni). However, in the usual 
regime of item [T] on page |4l they are always either real 
or purely imaginary because there is no set of quantum 
numbers for which parity-conserving and parity-violating 
contributions interfere with each other. 

To exploit the delta function in (l43l) we make use of 
the relation 

pldp, = Ei^Ef - M%dE, 

= (Mjv + T,UTy2MN +Ti)dT, (45) 


After integrating over T 2 we are left with the Ti integra¬ 
tion only. 


8 v-^ 

r^ = -E— / dT,p{T,,T2) 

A Jo 


K 1 K 2 J 


X (KiPitiK 2 P 2 ^ 2 d|A|KAAAfd) 


(48) 


where 


{KlpltlK2P2t2J\A\KAKNJ) 

= ■E((KipHlA2P2l2d|A|KAAArd) 

{K2P2hl^lPltlJ\^\l^Al^NJ))- (49) 

The direct matrix element is given by (1421) and the 
exchange one is obtained through the transposition 
(Ki,pi,ti) O- {K2,P2,t2)- 


B. Hypernuclei with open-shell cores and without 
recoil 


So far everything was done in the strict framework of 
relativistic physics. In what follows we will make use 
of analogies with nonrelativistic calculations. From pre¬ 
vious works [2ll - [^ done by our group, we know that 
to describe the hypernuclei with open-shell cores within 
the IPSM it is enough to do the following replacement in 
Eq.® 

^ ^ Ff (50) 

Ja 

where the spectroscopic factor is given by 

Jf 

y ((.riwciidjiwL (51) 

Jf 


As previously mentioned, to evaluate the spectroscopic 
amplitudes {JcWJjfjWJp), instead of employing the 
c.f.p.’s [13 that have been thoroughly used in, both non¬ 
relativistic [ 13 , and relativistic (4^48| calculations, we 
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use the second quantization formalism. In (1511) the sum¬ 
mation goes only over the values of Jp that fulfill the 
constraint \Jc — jN\ < Jf < Jc + Jn- The values for J/ 
and Jc are taken from experimental data and, for most 
hypernuclei of interest, are listed in Table I of Ref. [^ . 
The resulting factors Fj^ are listed in Table II of the 
same paper. 

Therefore, when the recoil effect is not taken into ac¬ 
count, the NMWD transition rate in open shell hypernu¬ 
clei reads 



jNt\t2 

K 1 K 2 J 



dTip{TuT2) 


X 


{KiPitiK2P2t2J\^\KAKNJ) 


.(52) 

T2=A,^-Ti 


We note that, while Eq. (H51) is only valid for doubly- 
closed-shell hypernuclei, Eq. (l5^ is valid for both closed- 
and open-shell hypernuclei. 


C. Inclusion of the recoil effect 

As seen above, when the recoil is neglected one can per¬ 
form first the full angular integration J dpi J dp 2 , lead¬ 
ing to a great simplification of the resulting expression. 
It is self-evident that this cannot be done anymore in the 
presence of the recoil energy 

Er = \Jy^\+p\+pl + 2piP2 cos 6 »i 2 , (53) 

where M/j is the relativistic rest mass of the recoiling 
nucleus. However, once the hypernucleus is unpolarized 
(and unaligned), there is no preferred axis along which 
to orient vectors. Therefore, we can choose to orient P 2 
with respect to pi and write 


Eor the sake of convenience we will work here with the 
nonrelativistic limit for the kinetic energy of recoil, i.e., 
with 


Tr = Er — M/j = 


pI + 2 pl + P 1 P 2 cos 6 >i 2 
2Mr 


Mat 


(Ti +T 2 - 277^008 012 ) (57) 


which we consider to be good enough for the present 
purposes. Moreover, we neglect the binding energy of 
the recoiling nucleus, and take = Miv(A — 2). The 
transition rate becomes then 


riv = ^ ^ / dTidT2dcosei2p{Ti,T2) 

K 1 K 2 J 

I _ 1 2 

X |(KiPitiK2P2t2T|A|KAKAfT)| 

X -Tr-Tx- T 2 ). (58) 


To perform the integration on T 2 we introduce an au¬ 
xiliary variable x, defined as T 2 = P 2 /= x'^, i.e., 

S{T2 + T,+TR-A,,)dT2 = ^ 

C-/2 m A 2 2a^COS^i2 /7^\ n T 

X d{x^ +Ti- - ^ VTT) ^xdx 

A —2 2xdx i-j., _,T 

= A-l\x+-x-\ ~ ^ ^ )]’ (59) 

where 


+ y/Ti cos 012 

X = - 


A-1 


(60) 


J dpi J dp2 ■ ■ ■ = J dpi J dpi2 ■ ■ ■ 

= J d4>i J d cos 01 J d(j)i 2 J dcos 0 i 2 --- ■ (54) 


Consequently, we can use (l2^ for integration on pi and 
(USD for integration on pi 2 , with the result that, as shown 
in Appendix [Cl instead of (l30l) we have now 


S'^ {piti,p2t2) 

= T. J 012 d(A,„ -T 1 -T 2 - Tr) 

K 1 K 2 J 

X |(piKitiP2K2t2d|A’"|/CAKjvJ)7 (55) 

Erom comparison with (l30l) one concludes that the results 
developed so far hold valid even when the recoil effect is 
included, as long as one makes the replacement: 


J p\dpipldp 2 5 {Aj^ - Ti - T2) • • ■ 

If 22 

-^2 dcos9i2Pidpip2dp2 

X d(A,„ -Ti-T2-Tr)--- . 


(56) 


Therefore 




3Ntlt2 

K\K2J 


X J dTidcos9i2dxxp{Ti,T2) 

X \{KiPitiK2P2t2J\A\KAKNJ)\ 

S{x — a;+) -I- d{x — x~) 

|a;+ — a;“| 

After integrating on x one gets 

jNtit2 

K 1 K 2 J 

X J dTidcos9i2^ p'(Ti,r 2 ,cos 0 i 2 , Ajjy) 

X \{KiPitiK2P2t2J\A\KAKNJ)\ 

}■ 


(61) 


L J X—^X~ 


(62) 

























where 

p'(Ti, T 2 , cos 6 * 12 , Ajjy) = 

_ x{A-2)p{T,,T2) _ 

v'Ti cos2 012 + (A - 2)(A - 1) - Ti{A - 2)2' ^ 

It might be useful to mention here that: 

• In the analogous nonrelativistic formulation, it has 
been shown numerically that the contribution cor¬ 
responding to the second term in Eq. (15^ is neg - 
ligibly small compared to that of the first 
Whether this also occurs in relativistic calculations 
must still be verified. 

• In the limit >1 —>■ 00 , the result (l52l) is recovered. 
Indeed, once 

- Ti, 

x~ becomes unphysical. Therefore, the only con¬ 
tribution comes from the first term in (1621) , and, as 
can be seen from ([63]), 

J dcosdi 2 p{Ti,T 2 , cos ^ 12 , Ajjy ) 

2p(Ti,r2). 

A. —> 00 

III. NUMERICAL RESULTS 

We present here our results for the NMWD rates of 
The recoil effect has been neglected since we have 
learned in our previous nonrelativistic calculations HH- 
[^. [ 2 ^ - 1 ^ that, although it is relevant for the energy dis¬ 
tribution of emitted particles in very light systems, such 
as s-shell hypernuclei, and is crucial for the angle distri¬ 
bution in general, it is less important for the integrated 
rates. Therefore Eq. (15^ has been used. 

Two approaches have been tested for the propagators 
A-^dx — y|), both based on the fact that the ranges 
of Yukawa-like baryon-baryon forces within hypernuclei 
depend not only on the intermediate meson mass but 
also on the baryon masses, as stated in [U, Appendix 
G], namely, 

RAl: This is the standard approach in nonrelativistic 
calculations [SElel, where the energy qo car¬ 
ried by the exchanged meson is constant and always 
smaller that the meson mass tom , having the value 
go = A/2 = 88.5 MeV. This implies that the fac¬ 
tor \/ — gp in (IT0|) is taking the place of the 
effective mass niM = ~ A^/4. 

RA2: This is the approach introduced in Refs. [13, H^, 
which is more appropriate for relativistic calcu¬ 
lations, where go is evaluated for each value of 
the kinetic energy Ti, with direct and exchange 
energies being respectively Qq = A + — Ti, 


and Qq = Ti — Once for the NMWD in 

A-hypernuclei the energy transfer is of the order 
50 — 150 MeV, go can be larger than ttItt and the 
factor can become complex. Therefore, 

in the case of the tt meson, besides making use of 
Eqs. ([T0|) and (ITT]) , one also needs Eq. (IT^ . We 
are particularly interested in this approach, since, 
as mentioned above, the transfer of energy in the 
NMWD of charmed nuclei can reach much higher 
values. 

Initial and final short range correlations (SRC) were 
included as in the nonrelativistic case, i.e., by making 
the substitution 

A^{r)^gf{r)A^{r)g,{r) (64) 

in the tensor expansion (1321) . where r = |x — y| and 



9f{r) = 1-joiqcr) (65) 

are, respectively, the initial and final SRC, with a = 
0.5 fm, P = 0.25 fm~^, and 7 = 1.28 fm, and gc = 
3.93 fm“^ [I 3 , m, m, I46l - l48j| . The dipole form-factor 
cutoffs Ajr = 1.3 GeV and Kk = 1.2 GeV are also the 
same as in these works. We present here two different 
sorts of comparisons involving our results for the decay 
rates of ^^C. First, in Table HIl the two relativistic cal¬ 
culations RAl and RA2 are compared with each other, 
and also with the analogous nonrelativistic (NR) calcu¬ 
lation using the RAl approach for the propagator. The 
NR calculation is analogous to the relativistic ones in the 
sense that it uses the same OME model, the same SRC, 
the same single-particle energies, and single-particle wave 
functions of a harmonic oscillator potential with size pa¬ 
rameter b = 1.60 fm, which gives the same root-mean- 
square radius for the initial hypernucleus. We show the 
decay rates r„ and Tp, the total one-nucleon-induced 
nonmesonic decay rates Tnm = r„ -|- Tp, and the ra¬ 
tios Tn/Tp within different OME models, namely, the tt 
and (tt, K) exchanges without and with SRC. Clearly, 
the relativistic calculations were evaluated in the labo¬ 
ratory frame of reference (LFR). Therefore, we confront 
them with NR calculations that also were done in the 
LFR. These, in turn, have been shown elsewhere [5^, [5^ 
to nicely agree with the NR evaluation within center- 
of-mass frame (CMF). It is not possible here to separate 
the decay rates r„ and r„ in the usual Block-Dalitz chan¬ 
nels 0 

a= iSo^^So, b= 3Po^iSo, c= 

d= ( 66 ) 

as one can always do in the CMF within the s-wave ap¬ 
proximation [24| . Therefore, we only show separate re¬ 
sults for the parity-conserving (PC) and parity-violating 
(PV) parts of the decay rates, which contain, respec¬ 
tively, the (a-l-c-f d) and (b-be-l-f) contributions. From 
Table El it can be concluded that: 
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TABLE II. Comparison between the nonrelativistic (NR) and 
relativistic results for r„, Lp, Tnm = r„ + Fp, and F^/Fp 
in for different OME models, i.e., the tt and n + K ex¬ 
changes without and with SRC. To allow a more detailed com¬ 
parison, parity-conserving (PC) and parity-violating (PV) 
parts of Fn and Fp are given separately. As explained in 
the text, two approximations were used for the propagators 
in the relativistic calculations, namely RAl and RA2. All re¬ 
sults are in units of the free A decay rate Fa = 2.50 x 
MeV. 


Model 

r(PC) 

-L n 

p(PV) 

A n 

pCPC) 

i p 

p(PV) 

i p 

Fn/Fp 

r 

A nm 

NR 

No SRC 

TT 

0.1084 

0.1592 

0.8717 

0.4017 

0.2102 

1.5410 

-k + K 

0.0286 

0.1851 

0.3748 

0.3550 

0.2929 

0.9434 

SRC 

TT 

0.0122 

0.1753 

0.8062 

0.4475 

01495 

1.4412 

n -\- K 

0.0161 

0.1945 

0.3822 

0.3985 

0.2697 

0.9913 

RAl 

No SRC 

TT 

0.1207 

0.1531 

0.6894 

0.3309 

0.2683 

1.2941 

■k + K 

0.0555 

0.2394 

0.3955 

0.3825 

0.3790 

1.0729 

SRC 

TT 

0.0969 

0.1226 

0.5210 

0.2447 

0.2866 

0.9853 

-k + K 

0.0563 

0.1648 

0.3503 

0.2691 

0.3569 

0.8405 

RA2 

No SRC 

TT 

0.1692 

0.2199 

0.7875 

0.4550 

0.3131 

1.6316 

-k + K 

0.0988 

0.3104 

0.4771 

0.5080 

0.4154 

1.3943 

SRC 

TT 

0.1454 

0.1865 

0.6209 

0.3619 

0.3377 

1.3147 

n + K 

0.1002 

0.2317 

0.4361 

0.3874 

0.4030 

1.1554 


1. Although it is clear that the inclusion of relativity 
sensibly affects the results, there is gross agreement 
between analogous nonrelativistic and relativistic 
calculations, both without, and with SRC. 

2. The SRC, while not crucial for some decay rates, 
can significantly reduce others, both in the nonrel¬ 
ativistic and in the relativistic cases. 

3. The decay rates r„ and Tp are both higher in RA2 
than in RAl. 

4. Relativity tends to make r„ become larger and Tp 

smaller, and this effect is more pronounced in the 
RA2 approach for the propagator. As a conse¬ 
quence the relativistic n/p ratio becomes signifi¬ 
cantly larger than the nonrelativistic one. There¬ 
fore, the relativistic approach, specially RA2, helps 
to solve the longstanding puzzle on the Tn/Tp ra¬ 
tio (See also Table Hill ) 

Then, in Table lllli are compared the present RA2 
calculations for with previous relativistic calcula¬ 
tions performed by Ramos et al. [d^, and by Conti et 


al. [T3, Idg for several OME models. The differences be¬ 
tween the three theoretical calculations are very large. 
We do not know the reason for such huge differences, al¬ 
though there are several possibilities. First, they could 
be due to the way in which the spectroscopic factors are 
evaluated, and how the states of the two emitted parti¬ 
cles are antisymmetrized and normalized. Secondly, dif¬ 
ferences can arise simply from numerical errors. To avoid 
those, we have checked step by step the entire relativistic 
calculation with the nonrelativistic one. 

In the same Table uni are shown the pertinent ex¬ 
perimental data produced by the KEK and FINUDA 
groups When confronted with the several rel¬ 

ativistic theoretical results, it is easy to discover that 
only the present evaluation within the t: + K + SRC 
model agrees well with the available data. In particu¬ 
lar, the good agreement for the ratio T^/Tp should be 
highlighted. The only significant discrepancy is with the 
experimental values for r„ and Tp obtained by KEK in 
Ref. [s^. However, there is agreement with the experi¬ 
mental value for Tp obtained by FINUDA in Ref. 571. As 
to the last column in Tables HIl and IIIIl it is important to 
remark that, while all the listed calculations include only 
one-nucleon-induced transitions, the experimental values 
include also eventual two-nucleon-induced contributions. 


IV. SUMMARY AND FINAL REMARKS 

Starting from the Fermi Golden Rule (|4]) we present 
in Section |n] a relativistic formalism to describe the non- 
mesonic weak decay of single-A hypernuclei within the 
framework of the IPSM, with the dynamics represented 
by the (tt, AT) OME model. First, in Subsection III Al we 
do this for hypernuclei whose cores have only closed sub¬ 
shells, and when the recoil effect is disregarded. Here, 
the Dirac plane waves are expanded in spherical partial- 
waves, the multipole expansion of the propagator is done, 
and the two-body matrix element is properly antisym¬ 
metrized with regard to the two outgoing nucleons. Ma¬ 
king use of the orthogonality condition (|2^ and exploring 
the energy conserving d-function, the six momentum- 
space integrals in Eq.(|3]) are reduced to one in Ea. (H51) . 
to be performed numerically. Next, the derived result is 
generalized to include hypernuclei with open-shell cores. 
This is done by means of the spectroscopic factors given 
by Eq. dsm, which are evaluated in second quantization, 
without recurring to the c.f.p. technique. In this way 
we arrive at Ea. (|5^ . Finally, in Subsection III Cl we dis¬ 
cuss the recoil effect, which is important not only for the 
evaluation of angular distributions of the pairs of emitted 
nucleons, but also for the study of single kinetic energy 
spectra in light and medium-weight hypernuclei. 

Numerical results for are presented in Section Ell 
Firstly, Table |ll] shows the comparison between analo¬ 
gous nonrelativistic and relativistic calculations of the 
transition rates r„ and Tp. The parity-conserving and 
parity-violating contributions are given separately. Such 
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TABLE III. Results for the nonmesonic decay rates r„, Fp and Vnm = r„ + Fp, and the ratio F„/Fp in for several OME 
models, namely, the tt and n + K exchanges without and with SRC. The present calculations with the approach RA2 for 
the p ropaga tors are compared with previous relativistic calculations performed by Refs, and with the experimental 

data All results are in units of Fq = 2.50 x 10“^^ MeV. 


Model 

F™ 

Tp 

Fn/Fp 

r 

nm 

TT 

Present (RA2) 

Ref. m 

Ref. [13, m 

0.39 

0.27 

0.89 

1.24 

1.32 

2.08 

0.31 

0.20 

0.41 

1.63 

1.62 

2.95 

TT-tSRC 

Present (RA2) 

0.33 

0.98 

0.34 

1.31 

Ref. [H 

0.06 

0.29 

0.20 

0.35 

Ref. [13, m 

1.80 

0.62 

0.34 

2.42 

n + K 

Present (RA2) 

0.41 

0.98 

0.41 

1.39 

Ref. [13, m 

1.24 

1.59 

0.78 

2.84 

TT -t K+SRC 

Present (RA2) 

0.33 

0.82 

0.40 

1.15 

Ref. [13] 

0.05 

0.36 

0.15 

0.41 

Ref. [13, El 

0.96 

1.42 

0.67 

2.38 

Experiment 

Ref. [531 

0.23 ±0.08 

0.45 ±0.10 



Ref. [13 

- 

- 

0.51 ±0.13 ±0.05 

- 

Ref. [13 

— 

— 

— 

0.828 ± 0.056 ± 0.066 

Ref. [13 

— 

— 

— 

0.953 ± 0.032 

Ref. [13 

- 

0.65 ±0.19 

- 

- 


a comparison is crucial so that one can rely on relativis¬ 
tic calculations. However, it had never been done before. 
The agreement is satisfactory, and the only difference 
worth mentioning is that the ratio Tn/Tp is apprecia¬ 
bly higher in the relativistic calculation, agreeing better 
with experiment than the nonrelativistic one, specially 
when the RA2 approach to the propagators is used. Sec¬ 
ondly, the present calculation is compared with two sim¬ 
ilar studies in Table Ell from where it is clear that the 
discrepancies are very large. Although there are impor¬ 
tant differences with our formalism, we could not find any 
justification to explain this. Done in the same table, the 
comparison of our calculation with the available experi¬ 
mental data is encouraging and favors the it + K SRC 
OME model. 

In short, we have achieved our goal of developing a 
reliable relativistic model for calculating the nonmesonic 
weak decay of A-hypernuclei, which can now be extended 
for similar weak decays in charmed nuclei. In parallel 
with this development, we intend to incorporate in our 
formalism the final state distortions of the outgoing nu¬ 
cleon waves induced by the interaction with the residual 
nucleus by making use of a relativistic optical potential. 
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Appendix A: Derivation of Eq. (1231) 

From the definition (fTOl) it follows that 

2 / ^dpj"^%'(ps|KTO)* (pslAt'm) • • • 

sm 

= 2j-^Sjj>{4TTf I dp^i'-''^(Z^is|jTO)y;*(p) 

sm (I 

■ ■ ■ 

= f dp^ l/;(p)lj4p).... (Al) 

















11 


Now, if we use the relation 

47r5]r,;(p)yi^(p)=P, (A2) 

A* 

we can solve the integral over the azimuthal angle to 
obtain 

2 J dpj~^(5jj/(ps|Km)*(ps|At'm) ■ ■ ■ 

sm 

= < 5 K„/( 47 r)^ y dcos0---. (A3) 

Appendix B: Relativistic single-particle wave 
functions 

The evaluation of the matrix elements of the NMWD 
is made in the context of the IPSM. This means that the 
A wave functions are those generated by spherically sym¬ 
metric mesonic mean fields. That is, in solving the Dirac 
equations for the single-particle level of A, one must use 
the meson mean fields from the nucleus. This is sim¬ 
ilar in spirit to the works of Ramos et al. [4^, |4^ , where 
single-particle bound-state wave functions are obtained 
by solving the Dirac equation with static Lorentz-scalar 
and -vector Woods-Saxon potentials. 

The radial bound-state wave functions F^ir) and 
Gi^{r) in (IMl) and corresponding energy eigenvalues £« 
for a single-particle state k for the N oi A are obtained 
by solving the following Dirac equations: 

( —I—'j Ff^ + (ek — D -I- S') G„ = 0, 

\dr r J 

^--)G^-(e.-V-S)F,= 0, (Bl) 

dr r J 

where the scalar potential S = S{r) is 

S(r) = M-f 5 „cr(r), (B2) 

with M = Mtv and for the N, and M = Ma and 

f/cr = g^ for the A; the vector potential V = V(r) for the 
nucleon is given by 

= 9^ ^o(r) +t^gpPo(r) + (t^ + l/2)eAo(r), (B3) 

with = 1/2 for the proton, = —1/2 for the neutron, 
and for the A it is given by 

V(r)=g:}uJo(r). (B4) 

The meson and Coulomb fields satisfy the following 

Klein-Gordon and Poisson equations 

(-V^ -f m^) cr = -g^ pf - 520-^ - gscr^ 
(-V^+m^JoJo=g^p^, 

(-V^ + m^)po = 1/2 gpps, 

-V^Aq = epp, 


with the densities given by 

K 

K 

K 

(B6) 

K 

where n// are the nucleon occupancies of the state k. 

The system of equations is solved by iteration follow¬ 
ing the scheme of Ref. [i^ : (i) we solve the Dirac equa¬ 
tions for given initial ansatze for the S and V potentials; 
(ii) the solutions for F{r) and G(r) are then used to solve 
the Klein-Gordon and Poisson equations and construct 
new potentials; and (iii) we put these into the Dirac equa¬ 
tions and cycle until convergence to a prescribed preci¬ 
sion is attained. Note that the nonlinear terms for the a 
field are put together with the scalar density p^ in the 
iteration procedure. 

The numerical values of the meson-nucleon parame¬ 
ters are those of the column NL3 of Table I in 
Ref. [13, and for the meson-lambda couplings are those 
from Ref. 0 (masses are given in MeV): 

g// = 10.2169, g// = 12.8675, gp = 8.9488, 
eV47r = 1/137, g^ = 0.464 g// = 0.481 

32 = -10.4307 fm-\ 33 = -28.8851, 
rua = 508.1941, rrii^ = 782.501, rrip = 763.000, 

Mat = 939, Ma = 1116.06. (B7) 


TABLE IV. Single-particle energies for and a^C. (Se e 
text.) Experimental values for are taken from Ref. 
and for a^C from Ref. [tII]. All values are in MeV. 



Calculated 

Experiment 

P lSl/2 

-38.53 

-34 

P lP3/2 

-13.52 

-15.96 

n lsi /2 

-42.03 

-37 

n lp3/2 

-16.65 

-18.72 

A lSi/2 

-11.59 

-10.79 


In Table IVl we present the single-particle energies 
for ^^G and a^C. Note that these results are obtained 
without adjusting any parameters to fit experimental 
numbers. Clearly, a reasonable description of the exper¬ 
imental single-particle energies is achieved. Of course, a 
better description could be obtained by fine tuning the 
parameters, but for the purposes of the present paper 
such a refinement is not necessary. 


(B5) 
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Appendix C: Derivation of Eq. (1551) 


to obtain 


Here we demonstrate the result (1551) starting from the 
definition (H?]) for S'^{piti,p 2 t 2 ), i-e., 

S^{piti,p2t2) = ^ f dpidp2S{Ajj^ - Ti - T 2 - Tn) 

m\mN 

S 1 S 2 

X |7W’"(piP2SiS2tH2jjVTOjvjATOA)P- (Cl) 

Using the expansion (I25|) and making the change of vari¬ 
able P 2 —!> P 12 as explained in Eq. (IMl) . we are free to 
perform the pi integration according to Eq. (|2^ . and are 
left with 


S'^{piti,P2t2) = 

(4^)2 ^ J #i25(A,^ _ Ti - T 2 - Tr) 

m\mN 

K-im\S2 

X I ^ (Pl2S2|«:2W2) 

K2m2 


X {piKimitip2n2m2t2\A'^\KKmRKN'mN) 


2 


(C2) 


Then we do angular momentum algebra as in (I29L 


S'^ [piti,p2t2) = 

{Airf Y. f dpi2S{A,^ - Ti - T 2 - Tr) 

S2 

X Y (Pl2S2|K2m2)*(pi2S2|K2m2) 

K2m.2 

X Y (il”^U2W2| JM)(jlmlj2^7^2| JM) 

J 

m\M 

x{piKxtiP2K2t2J\A'^ \rrKn jy 
x(piKlUP2K2f2^|A’"|KAKAr J). (C4) 

Due to the relation 

Y ^nmij2m2\JM){ji'mij'^m'^\JM) = 

mi M 
J2 

’ (C5) 

J 2 

this reduces to 


S^{piti,P2t2) = 

(4^)2 ^ I dpi 2 ^(A,„ -T 1 -T 2 - Tr) 

m^rriN 

KimiS2 

X Y, (Pl2S2|K2W2)*(piKitiP2K2<2-^|A’"|KAKArJ)* 
K2m2 
JM 

^{jimij2m2\JM){jRmRjMmM\JM) 

X Yj (Pl2S2|K2m2)(piKitiP2K2i2^'|A’"|KAKAr^') 

J'M' 

x(jimij2W2| J'M')(jATOAjArTOAr| J'M') , (C3) 


S^{piti,P2t2) = 

(47r)2 ^ / dpi2^(A,, - Ti - T 2 - T«) 

K2m2 -^2 
K. 2 S 2 

x(pi252|/C2m2)*(pi2S2|«^W2) 

X ^ J^(piKitiP2K2f2>/|A’^|KAKv J)* 

X (piKitiP2K2^2^|A’^|KAKAf^)- (C6) 

Finally, using (l23l) one gets (l5^ . 
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